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Exploring qubit-qubit entanglement mediated by one-dimensional plasmonic 

nanowaveguides . 
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We exploit the qubit-qubit coupling induced by plasmon-polariton modes in a one-dimensional 
nanowaveguide to obtain various two-qubit entanglement situations. Firstly, we observe three phe- 
nomena occurring when preparing the initial state of the system and leaving it freely relax: spon- 
taneous formation of entanglement, sudden birth and revival. Then, we show that plugging a laser 
to each of the qubit, the system arrives to an steady state which, depending on their inter-qubit 
distance, can also be entangled. For this situation, we also characterize the quantum state of the 
system showing the entanglement-purity- diagram typical for two-qubit systems. 

I. INTRODUCTION 



Entanglement is one of the most striking features of 
quantum mechanics. From the fundamental point of 
view, the non-separability of systems in product states 
of their components is a very interesting phenomenon as 
it has no classical analogue. From the practical point of 
view, entanglement is a key element of quantum cryp- 
tography, quantum teleportation, and other two-qubit 
quantum operational. Although it was first studied and 
exploited for atomic-like systems, due to the advances in 
semiconductor physics and their promising scalability in 
fabrication important progress have been made in this 
direction and short-range entanglement between spin or 
charge degrees of freedom in quantum dots, (QDs), nan- 
otubes, or moleculea^H^ has already been observed. How- 
ever, in order to get some long-distance entanglement 
one need some kind of interaction mediating between 
the qubits is needed. Usually photons^^ have played 
this role, although there have also been some other pro- 
posals, for instance using NV centers in diamond^ or 
more recent theoretical proposals to use one-dimensional 
plasmon-polariton (hereafter simply labeled as plasmons) 
nanowaveguides to achieve steady-state entangled two- 
qubit states^SEH. 

Surface plasmon-polaritons(SPP) are confined 
electromagnetic modes occurring in metal-dielectric 
interfaces'^. In the case of channel plasmonic waveg- 
uides(PWs) like the one of Fig. fl] which is the case 
that we will be concerned with here, these modes have 
one-dimensional behaviour and have been proposed as 
useful elements for the future generation of photonic 
circuits^. The interest in using them for quantum op- 
tics applications'^ rose considerably when the coupling 
between a quantum emitter and a single plasmon mode 
was observed'^. These works show that the /3- factor, 
which measures the fraction of the emitted radiation 
that is captured by the propagating mode, can be close 
to 1 due to the subwavelength nature of the plasmon 
field. Extension of this formalism to the case of two 
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FIG. 1. (Color online) (a) Two qubits interacting with a 
plasmonic waveguide, in this case a channel waveguide, (b) 
Scheme of levels, couplings, and decays in the particular case 
where uji =^ 0)2 = (jOq and 711 = 722 = 7- 



emitters placed close to this PWs it was also showrP^ 
that one could get a modulable coupling between them, 
as a function of the inter-qubit distance, which could be 
used for generating entanglement between them''. 

The paper is organised as follows: in section |TI| we 
describe the theoretical framework that we will use for 
our calculations, i.e. the master equation formalism; in 
section IIIII we show the results where we consider that 
our system is initialized at some given state and then let 
it freely relax. Finally we consider a situation where the 
qubits are coherently pumped in section IV , and conclude 
afterwards in section IVl 



II. THEORETICAL FRAMEWORK. 

From the quantum optical point of view, we are going 
to see our system as two qubits, which are interacting 



with the plasmonic field in the PW, characterized by the 
continuum of bosonic modes. When the coupling be- 
tween the qubits and the continuum of modes is weak, 
then one can safely perform a Born-Markov approxima- 
tion, and trace out over the degrees of freedom of the 
SPP. This approximation leads to he dynamics of the 
density matrix p for two qubits described by a master 
equation^^-^ as follows: 



dtp = i[p,H]^ Y, ^-f(2cj,pcj]-cj\cjjp-pcj\cj,) (1) 

where crj,cr^ are the raising and lowering operators for 
each qubit, and where H^ taking two qubits with charac- 
teristic frequency uOi is given by: 

H = hY {uJiajai +^i2(cr|cr2 +crjcri)). (2) 

i=l,2 

The SPP play a double role: in the coherent part of the 
dynamics an effective interaction between the two quits, 
also present in cavity QED, is provided by the exchange 
of virtual plasmons: 
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^/m[/i^G(cj,ri,r2)/i2] 



duj, (3) 



where G(ri,r2) is the off-diagonal Green's function 
describing the electromagnetic interaction between two 
dipole moments /ii^2 of frequency uq. The rates of diag- 
onal and off-diagonal dissipations are given by 



2^2 
^ij = — |-/m[/i*G(a;o,ri,r^)/i 



J J' 



(4) 



with z, j = 1, 2, 7i2 = 7 21 an d we will label ju = 7^. 

It was recently founcP^EU that when the propagating 
plasmon supported by the PW is the dominant channel 
for emission (i.e., large /3 factor), a very good approx- 
imation for the total Green's function can be obtained 
by only including its plasmon contribution, G(ri,r2) ~ 
Gp(ri,r2). In this way, analytical expressions for both 
gi2 and 712 can be easily derived: 

7. 
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= 7/3e-^/(^^)cos(fc^,d). 



(5) 



where d is the separation between the two qubits, kp 
and L being the wave-number and propagation length 
of the plasmon, respectively. The crucial point of Eq. 
[5] is that the 7r/2 phase shift between the coherent and 
incoherent parts of the coupling mediated by plasmons 
enables to switch off one of the two contributions while 
maximizing the other by just tailoring the inter-qubit 
distance, opening the possibility of modulating the degree 
of entanglement as we will see in the next sections. 



III. SPONTANEOUS DECAY. 

In order to solve the dynamics given by the mas- 
ter equatior^T] we need to project in some basis this 
equation and then solve the set of differential equa- 
tions for different initial conditions depending on our 
initial state. A possible basis to represent the dy- 
namics (IT]) is the one given by {|0) = |^i,^2),|l) = 
|ei,^2),|2) = |^i,e2),|3) = |ei,e2)}, where gi/ci la- 
bels the ground/excited state of the i-qubit. However, 
for a better understanding of the physics behind entan- 
glement it is convenient to change to the basis {|0) = 
1^1, ^2), |±) = -^{\l) ± |2)),|3) = |ei,e2)}, where two 
maximally entangled states |±) appear explicitely. This 
basis also allows to understand more clearly the dynam- 
ics of the system as diagonalize the coherent part from 
the dynamics. In figure [l](b) could be observed that de- 
pending on both the sign and absolute value of the cross 
decay term 712 one of the states |±) could be decoupled 
from the dynamics of the rest of states. 

Once the density matrix p{t) is obtained by numeri- 
cally solving Eq. (IT]) represented in the mentioned basis, 
the entanglement of the two qubits is quantified by means 
of the concurrence^^, C = [max{0, v^Ai — >/A2 — a/As — 
V^}], where {Ai, A2, A3, A4} are the eigenvalues in de- 
creasing order of the matrix pTp*T, with T being the 
anti diagonal matrix with elements { — 1,1,1,— 1}. The 
two main ingredients of our problem affect C in a differ- 
ent way. The coherent coupling gi2 produces oscillations, 
whereas the effect from the cross-decay 712 is much more 
important non-oscillatory contribution to C. 

Now we will describe two situations where the system 
is initialized in a given state from which it decays sponta- 
neously. In the first one, the system is initially prepared 
in a one-excitation unentangled state and in the second 
one a doubly-excited state. 



A. Spontaneous formation of entanglement. 

In order to analyze the spontaneous formation of en- 
tanglement in the system, one can initially prepare the 
system, for instance in the state |1), by a 7r-pulse on one 
of the qubits. Then the only non-zero elements of the 

initial density matrix are: p++(0) = p (0) = p+_(0) = 

p ^(0) = 1/2, so the dynamics is confined to the sub- 
space spanned by these three vectors: {|0) , |±)}. Apart 
from the populations, the only non-zero elements in the 
dynamics are p±^(t). The expressions for the non-zero 
density matrix elements in this case are given by: 



P++{t) = -7+(p++W) 

P— W = -7-(P— W) 

-+{t) = -{j-2igi2)p-+{t) 



(6) 



_ 71+72 

takes the form: 



where 7± = ^^^^^ ± 7i2- In this case, the concurrence 



C{t)^/[^^^^ 



_(t)]2+49[p+_(t)]2 



(7) 



and solving Eq. |6] we arrive to the fohowing simple 
expression for concurrence: 



C{t) 



-It 



sinh^(7i2t) +shi{2gi2t) 



(8) 



where we can see there are two different contributions: 
the one from the coherent coupling gi2 which is respon- 
sible for the oscillations, and a non-oscillatory contribu- 
tion corresponding to 712 which produces an asymmetric 
population of the |±) states responsible for the lengthen- 
ing of the lifetime of the entanglement in this particular 
situation. The dynamics of C and its dependence with 
the inter-qubit distance d is shown in figure |2] for a PW 
with /3e~^/^^^^ = 0.94. For all d one could observe an 
fast spontaneous formation of entanglement, however for 
d = nAp/2 with n = 1, 2, ... its decay is much slower. This 
is mainly due to the large asymmetry in the timescale of 
the two cascades in Fig. [T] which produces a large imbal- 
ance of populations and therefore creating entanglement 
as predicted in Eq. [7[ 



B. Sudden birth and revival of entanglement. 

Since one of the goals of this paper is to show all the en- 
tanglement phenomenology of our particular system we 
are going to focus now in a different situation where our 
system is initially prepared in a non-maximally entangled 
state: \^q) = y^ |0) + ^JT^ |3), where < g < 1. Con- 
sequently the initial non-zero elements of the density ma- 
trix are: poo(O) = q, ^33(0) = 1 - ^, Po3 = \/<7(l -q)- So, 
as in the previous case, the rest of the density matrix will 
remain zero except for the symmetric and antisymmetric 
populations that will build up during the evolution. The 
expressions for the non-zero elements are given by: 



P^z{t) = qe-'^' 



7- 
7+ 



oi-t 



o-r+t _ I 



Pm 



(t) = ^q{l - g)e-(7-2ic.o)t 



(9) 



After a tedious but simple algebra you arrive to 
the following expression for the concurrence C{t) = 
max{0,Ci(t),C2(t)}, where: 

Ci(t) = 2|po3(t)|-|p++(i) + P-(t)| 
C^it) = -2^poo{t)p33{t) + \p++{t) - p—{t)\ (10) 

This double criterion for entanglement leads to a very 
interesting and well-known phenomena which has been 
called as revival of entanglement^^. In figure [s] we have 
plotted the concurrence for a system with q = 0.5, which 
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FIG. 2. (Color online) Panel (a) Density plot of concurrence 
as a function of time and distance between qubits when the 
system is initially prepared at the state |1). It corresponds 
to two qubits entangled by means of a SPP waveguide with 
^g-d/(2L) ^ Q_g4^ ^^^Yi Xp = 640nm and 71 = 72 = 7.The 
spontaneous formation of entanglement is observed. Panel 
(b) show the concurrence dynamics for d = Xpi in black and 
the |+)/|— ) population in dotted blue/dashed red line re- 
spectively. 



is initially entangled. It could be seen that after a short 
time, of the order t/j = 4 entanglement disappears for 
a time, but latter it, counter intuitively, revives. The en- 
tanglement revival is clearly related with 712 as we could 
see that in the distances where it is small, the revival 
does not take place. Counterintuitively, as the {|1) , |2)} 
(or |±)) are initially uncoupled, and remain uncoupled 
forever, the coherent coupling gi2 of the dipole does not 
play a significant role in the dynamics of the system, as it 
can only change one excitation, but not two. Quantita- 
tively the revival is very subtle as it only gets a maximum 
entanglement of around 3%. 



Considering the limiting case where g = 1, so that our 
initial state is given by |$o) = |3), which is unentangled, 
then the only non-zero element of the initial density ma- 
trix is given by P33 = 1. Therefore the p{t) remains 
diagonal during its evolution, being the populations the 



only non-zero elements: 

p-{t) 



P33{t) = e 
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(11) 

One would expect then that no entanglement could 
build up from this initial state, however, after some sim- 
ple algebra you arrive to the following expression for the 
concurrence C{t) = max{0,C2(t)} where: 



C2(t) = |p++(t) - p__(t)| - 2v/pooWP33W (12) 

So if the difference of population between the symmet- 
ric and antisymmetric states overcomes the second term 
in Eq. [T2J some concurrence may build up. Again, the 
crucial element is the collective decay rate 712 and not 
the coherent coupling ^12, which does not play a role. In 
figure |4] this phenomenon could be observed for a given 
set of parameters (cf. caption): no concurrence occurs 
for the initial times and suddenly a 2% is built. This 
phenomenon usually called sudden birth'^^ ^^ of entangle- 
ment, is again subtle in our system but still observable. 



IV. STEADY-STATE ENTANGLEMENT. 

Up to now, we have seen that it is possible to observe 
finite-time entanglement in different situations. Never- 
theless, one is usually interested in having a stationary 
state with some degree of entanglement. In order to 
achieve a stationary-state a continuous pumping is re- 
quired. The pumping of the qubits may have incoher- 
ent nature^ or one could use coherent excitation, i.e. 
a laser field. We are going to assume a situation for 
sufficiently separated qubits, where the stationary state 
can be modulated by acting independently on each qubit 
with a laser of Rabi frequency Qi. Therefore a new term, 
"^^ hQi{al-\-ai) must be included in the coherent dynam- 
ics, i.e. in Eq.Q. 

As we recently showed^ one could find steady-state 
entangled situations playing with the laser intensities. In 
Fig. [5J we have plotted the steady-state concurrence for 
the typical set of parameters that we have been consid- 
ering(cf. caption) within the paper as a function of the 
inter-qubit distance and for different pumping configu- 
rations. In this figure, the pumping on the first qubit 
is fixed: Qi = O.I7, while the second laser intensity Q2 
ranges from —Qi to l^i, going then from an antisym- 
metric pumping with Qi = —Q2 where one can observe 
that peaks in the concurrence appear for d close to an 
even multiples of Xpi/2. As expected, in the symmetric 
pumping situation (l^i = ^^2) the peaks appear for odd 
multiples of Xpi/2. In the intermediate region of asym- 
metric pumping, Qi ^ ^2^ the periodicity is broken and 
both even and odd multiples of Xpi/2 appear. 
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FIG. 3. (Color online) Panel (a) Density plot of concurrence 
as a function of time and distance between qubits when the 
system is initially prepared at the state |$o) defined as in the 
text with q = 1/2. It corresponds to two qubits entangled 
by means of a SPP waveguide with /3e~'^^^'^^^ = 0.94, with 
Xp = 640nm and 71 =72 =7. The so-called revival of entan- 
glement is observed. Panel (b) show the concurrence dynam- 
ics for d = Xpi in black and the |+) / |— ) / |3) population in 
dotted blue/dashed red/long-dashed green line respectively. 



This possibility of modulating the steady-state of a 
system by changing the pumping and assisted by dissi- 
pation is connected with the ideas of dissipative com- 
putation and state-engineering recently proposed^^ and 
experimentally demonstrated^. 



A. Entanglement-purity diagram. 

Up to know we have only worried about the quantifi- 
cation of the degree of entanglement using the Wooters 
concurrence. However, as we are dealing with a dissi- 
pative environment, the quantum state of the system is 
never pure and its degree of purity is very relevant for its 
possible application in quantum information protocols^ 
As we have already seen, we use a density matrix de- 
scription which is a very powerful formalism to describe 
mixed states. One standard measurement for character- 
izing the degree of mixture of a system with a given p 
density matrix is the linear entropy'^^ defined by: 



5i = -(l-IV(p2)) 



(13) 
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FIG. 4. (Color online)Panel (a): Density plot of concurrence 
as a function of time and distance between qubits when the 
system is initially prepared in the state |3). It corresponds 
to two qubits entangled by means of a SPP waveguide with 
^g-d/(2L) ^ Qg^^ ^.^j^ ^^ ^ 640nm and 71 = 72 = 7. The 

so-called sudden-birth of entanglement is observed. Panel 
(b) show the concurrence dynamics for d = Xpi in black and 
the 1+) / |— ) / |3) population in dotted blue/dashed red/long- 
dashed green line respectively 
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FIG. 6. (Color online) Concurrence-Linear entropy diagram 
for a system with parameters /S = 0.94, L = 2/im with each 
point corresponding to different inter-qubit distance ranging 



from d/Xs 



0.5 — 1.5, whereas the three different color 



^spp 



correspond to three different pumping configurations: Qi 
0.157 and Q2 = (red), Qi = Q2 = O.I7 (blue), and Qi 
— Q2 = O.I7 (green) 



When one has a pure state then Tr(p^) = 1 so the 
linear entropy becomes Sl = whereas for a maximally 
mixed state Tr(p^) = 1/4 so that 5'l = 1. However, the 
most interesting regime for us is what happens in between 
when some entanglement is present but the purity of your 
system is not perfect. It is well-known that if the den- 
sity matrix describing a two-qubit system with a certain 
degree of mixture, then there is a maximum degree of en- 
tanglement this system may achieve by means of unitary 
transformations. These states are usually called Max- 
imally Entangled Mixed States (MEMS) and were first 
proposed by Ishizaka and Hiroshima^^ and they occupy 
a region the concurrence-linear entropy diagram^^ shown 
in black in Fig. [6] We have also included the points re- 
sulting from the calculations in our system for the three 
different pumping configuration that we aforementioned 
in the previous section (see the caption for details) so that 
one could see how far from the optimal entanglement is 
our state for a certain degree of purity. 



CONCLUSIONS. 



FIG. 5. Steady state concurrence as a function of the sepa- 
ration between two equal qubits for /3 = 0.94, L = 2/im and 
for a pumping configuration with Qi = O.I7 and Q2 ranging 
from — Qi to Qi. Black regions correspond to higher values of 
concurrence while the white regions correspond to zero con- 
currence. 



We have used the properties of the qubit-qubit cou- 
pling induced by plasmon-polariton nano waveguides ^^Hl 
to explore a wide range of entanglement situations: start- 
ing from unentangled states we have shown the possi- 
bility of spontaneous formation, sudden birth, and re- 
vival of entanglement. Using a pumped configuration 
scheme, we have been able to achieve steady-state mod- 



ulable entanglement and characterize the quantum state 
in the concurrence-entropy diagram, comparing them to 
the MEMS. The possibihty of modulating entanglement 
open the possibility of implementing the ideas of dissi- 
pative quantum computation^ in this type of plasmonic 
systems. 
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